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loss predictions are produced by computer methods. An attempt 
is made to define the standard atmospheric conditions along 
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normal and super-refracti ve propagation conditions. 
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I, INTRODUCTION 



Numerous methods of estimating tropospheric propagation 
path loss are currently found in the literature. These 
include both graphical and computer techniques, of which the 
graphical methods are most widely presented. In conducting 
this particular study a program dealing with tropospheric 
path loss prediction (TROPOPLOT) was chosen as a computer 
method of analysis. This program was originally published by 
ESSA in a technical report [Ref. 6] and was subsequently 
modified for use at the Naval Postgraduate School [Ref. 7]. 
Tests were conducted by Longley and Reasoner [Ref. 8], in 
which the computer-predicted results of the program were 
compared with empirical data taken over a number of different 
paths. The results showed that TROPOPLOT can provide 
reasonably accurate predictions of path loss within certain 
constraints. 

The paths selected for analysis in this study are located 
in the Puget Sound, Washington area. They comprise both 
existing and proposed communications links connecting the 
underwater range facilities operated by the Naval Torpedo 
Station at Keyport, Washington. The range operation centers 
are currently located at Winchelsea Island, British Columbia 
and on Zelatched Point near Dabob Bay. The present communica- 
tions requirements include on-range communications with both 
range control vessels and the submarines and aircraft conduc- 
ting tests, the capability to monitor telemetry information. 
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and the maintenance of telephone links with other sites. 
Proposals for expansion of the range facilities as the 
TRIDENT program progresses include possible at-sea range 
facilities and digital data transmission to a central computer 
processing center located at Bangor or Keyport. These pro- 
posals would increase the bandwidth requirements of the links 
and thus necessitate a change from the present VHP frequency 
range to an S-band (2 GHz) line of sight lijik. 

The terminal and repeater sites are shown in Figure 1. 

The existing VHP link consists of paths from Winchelsea to 
Lookout Mt. (repeater) to Bangor and from Zelatched Point to 
Bangor. The proposed link at 2 GHz runs from Makah to Striped 
Peak (repeater) to Mt. Constitution (repeater) to Gold Mt. 
(repeater) to Bangor or Keyport, as well as from Zelatched 
Point to Bangor or Keyport. Each of these paths is computer- 
analyzed in Section II. 

In addition to computing path loss, the historical 
meteorological data for this area was examined in an attempt 
to characterize the "typical" propagation conditions found in 
this vicinity. This was accomplished by means of the 
Integrated Refractive Effects Prediction System (IREPS) which 
provided output based on data taken from upper air sounding 
stations in the Puget Sound area. This information could also 
be useful in the determination of link service probability. 
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Figure 1 

NAVTORPSTA Keyport Ranges and Vicinity 
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A listing of a computer program designed to plot path 
profiles and to aid in the determination of several input 
variables is included in the appendices, as is an explanation 
of the modifications required in order to use detailed terrain 
profile information in the TROPOPLOT program. 
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II, COMPUTER-AIDED PATH LOSS PREDICTION^ 

A, PROGRAM DESCRIPTION 

TROPOPLOT provides as an output a measure of the long- 
term attenuation over a given path. The program was primarily 
designed for use in cases where detailed terrain information 
was not available but can be used accurately with user- 
supplied data after a slight modification. The required 
input parameters include frequency, antenna heights, antenna 
gains, line losses, receiver sensitivity, path length, 
polarization, transmitter power out, surface ref racti vi ty , 
conductivity, permittivity, and a measure of the terrain 
roughness Ah. These parameters are for the most part well 
defined in Refs. 6 and 7 and are easily determined except in 
the case of user-supplied terrain profile information. For 
this sort of input there is some ambiguity concerning the 
selection of values for a, e, Ah, and the effective antenna 
heights. 

In considering a path which includes both poor ground 
surfaces (e=4, a=.001) and sea water surface (e=81, a=5) the 
particular value to use as input is not specified, although 
a completely oversea path requires an adjustment to the 
program not available in the version used for this study. 

Thus the values of s and a chosen were based on the corres- 
ponding value for that portion of the terrain that constituted 
the largest portion of the dominant reflecting plane between 
the transmitter and receiver. 
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The path geometry considered is shown in Figure 2 and 

indicates h , « as the structural transmitter (receiver) 
gl ,2 

height above ground. This input parameter is converted into 
an effective height, h^.| 2 the program. It should be 
noted that Longley and Rice consider two cases: one for 
random antenna siting in which structural heights, hg.j ^ 
considered equal to the effective heights, hg.j 2 > and another 
in which the sites are carefully selected, as in the case of 
radio relay links, and the effective heights are larger than 
the structural heights as shown below: 

hgi ,2 = ^gi,2 kexp(-2hg^ ^ 2 /^h) meters 

The variable k is considered to have a maximum value of 
50 as determined by the author's study of varied terrain 
conditions. The method presented below for the computation 
of k is valid only for antenna heights less than or equal to 
10 meters. In Callaghan's version [Ref. 7] this determination 
was : 



k = 1 + 4sin(nh , o/lO) for 0 < h 1,2 < 5 

gl »2 ' - g - 

= 5 otherwise 

Since most of the antennas in the links under study had 
structural heights in the neighborhood of 50 meters, it is 
obvious that predictions based on strictly hg-j 2 would in this 
case, produce an erroneous (higher than normal) prediction of 
path loss. This factor is recognized by the authors as an 
area requiring better definition and as a primary source of 
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Figure 2 
Path Geometry 
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prediction error. For this reason the modifications to 
TROPOPLOT which are detailed in Appendix B were necessary. 

The effect of varied antenna heights in the input of the 
unmodified program is shown in the output data presented 
later in this section. 

The description of the terrain is accomplished by a 
statistical quantity, Ah, which is defined in Ref. 6 as "the 
asymptotic value of the interdecile range, Ah(d), of terrain 
heights above and below a straight line fitted to elevations 
above sea level. The parameter Ah(d) is calculated at fixed 
distances and its median value usually increases with path 
length to Ah." For a single path profile these definitions 
are not adequate. The asymptotic value. Ah, used in this 
study as determined by taking the interdecile range of the 
difference between a straight line fitted to the path profile 
points and the path profile points themselves. The results 
of this estimation of Ah correlate favorably with those 
values listed in Table 1 [Ref. 7] which shows estimated values 
of Ah for particular types of terrain. Figure 3 shows an 
example of the straight line fit to the terrain profile between 
Lookout Mt. and Bangor. The details of the computation are 
found in Appendix B. 

Having thus considered the variations in the input 
parameters the program constraints should be noted. TROPOPLOT 
is designed for use within the following constraints; 
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Input Values for Various TROPOPLOT Runs 
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Figure 3 

Examples of Least Square Line Fit to Terrain 
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Parameter 



Range 



frequency (f) 

antenna height (h^^ 2 ) 

distance (dist) 

surface refractivity (N^) 



20 to 40,000 MHz 
0.5 to 3000 m 
1 to 2000 km 
250 to 400 N-units 



The antenna siting is subject to the following conditions; 

1. The angle of elevation, 8^^ 2 » O'*" each horizon 

ray from the horizontal should not exceed 12°. 

2. The distance from each antenna to its horizon 

(d^l 2 ) should not be less than 1/10 or more 
than 3 times the corresponding smooth earth 
distance 

These limitations can be ignored if the values of d, , 

h ,2 

and 2 entered directly into the program. As is the 

case with Ah and h^^ 2 method of accomplishing this is 

contained in Appendix B. 

Having resolved, or at least indentified, the potential 
ambiguities in the input parameters, it is useful to briefly 
describe the output format of the program. As shown in 
Figs. 4-7, the output of TROPOPLOT consists of two tables 
containing transmission loss versus distance and signal 
strength versus distance, a plot of the transmission loss 
information and a printout of the calculated and input values 
of many of the program parameters. Some of the quantities in 
this figure may require some explanation as to their meaning 
or derivation. These include the variables AE, AES, AED, MS, 
^ MD, ADX, DX, K1 , K2 , and ALS. 



22 



« 



'Ifir 



k .» t<9 

m 

I <103 

I 

• I 
I 

I 

>1 e i ' 

'9 ? 



»f 



FREOUENCY OF SIGNAL - F- INE6-HERT2) 100.00000 



o) 



o 


o 


o 


o 


o 


o 


.J— 


in 


OO 




>r 


e- 


1 ^ 


• 


• 


• 




• 


OOCJ 


nT 


>r 


m 


in 


m 


Zw 




r* 


p>^ 


e- 


f'- 


< 












of 


I 


1 


1 


1 


1 


»— 


1 

1 


1 

1 


1 

I 


1 

1 


1 

1 


UJ 


1 

1 


1 

1 


1 

1 


1 

I 


1 

1 


o 


1 


1 


1 


1 


1 



z-^ 


O 


o 


o 


o 


o 


<z 


O 


OJ 


o 


o 


o 




O 


o 


o 


o 


o 


1/1— » 


CO 


gj 


vT 


eg 


o 




* 


• 


• 


• 


• 


Q 




m 


g3 




00 




o 


gj 


gj 


gj 


gj 



O O ^ (J' 00 o 

O ft) m CO ^ nT 



••••••• 

o <o o o r- 

r» p«. f%- r«- p» r- 




o o o o o o o 

o o o o o o o 

0 0 ^ 0 0 0 0 

CO o %r eg O 00 o 

• ••*••• 

00 O' o ^ <M fT^ 

•O o r*. 



inNj-m.HOo>r-in 

• ••••••* 

COCOCOoOV'O'CT'O 

r^p^r*fs»/v,r».pH,'D 




OOOOOOOO 

*^0000000 

OOOOOOOO 

'3*<\iOcogjvrego 

'TIAO>0^-C00'0 

f^r^r~r*r^e-f^ao 



uo 

uo 

0 
— 

1 o 

t/>Q 
Z — 
< 

QC 



(\| O O' 

— I <r r- 



00 CO 

o <J 



CO O' O' 

<J o o 



00 ^ o CO 

^ nt fv. o 

• • • • « 

O' O o o — • 

o e- f'* e- 



o eg m nT o 

vj' e- o vO 

• • • • • 

^ ej pvj {\j 

e- p^ p^ f'- 



P^ CO O' 

O' eg in 

• • • 

eg m m 

r- r-. 



O' O 
00 eg 






o 

Z-' 

<2. 

to*— 

O 



r 



o 


CJ 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


00 


o 


>r 


eg 


o 


CO 


o 


• 


• 


• 


• 


• 


• 


• 


00 


O' 


o 


•■g 


eg 


ej 


re 


>r 


vT 


in 


in 


m 


in 


in 



o 


o 


1 

o 


1 

i-j 


1 

o 


1 

‘ j 


1 

o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


'T 


eg 


o 


X 


o 


nT 


eg 


• 


• 


• 


• 


• 


t 


• 


vT 


m 




X 


e- 


X 


CJ' 


m 


in 


in 


in 




in 


m 



O O o o o o 

O O O O O o 

O O o o o o 

O CO o vl“ eg o 

••••*• 

OO^egrr^vT 
O O go O O 



I 

I 

I 

< 

H- 

< 

O 



< 

QC 

Oi 

(U 



O’ 

UJ 

00 

o 



I 

I 



oo 

oo 



O 


m 




o 


m 


sT 


n 


-t 




r- 


vj* 




e- 


re 


O' 


m 


O 


in 


O' 


vj- 


X 


1 X 






m 


O' 


re 






<r 


X 


fSJ 


gj 


O' 


re 


gj 


o 


>r 




O 






OOQ 


o 






-g 


eg 


eg 


re 


cn 


re 


>3“ 


>3“ 


'T 


in 


in 


o 


gj 


X 




e- 




Z — ’ 


gj 


gj 


gj 


gj 


gj 


gj 


o 


>«j 


gj 


O 


gj 


gj 


gj 


gj 


X 


gj 


gJ 


gj 


gj 


O 


or 


1 




i—t 


i—t 














*-• 


-g 




<-g 














J— 


1 

I 

1 

1 


1 

1 

1 


1 

1 

1 


1 

f 

1 

1 


1 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

I 

F 


1 

1 

< 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

I 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


Y 

1 

1 


UJ 

O 


f 

1 

1 


1 

1 

1 


I 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

i 

1 


1 

1 


1 

1 

1 


1 

1 


1 

1 


1 

1 

1 


1 

1 


I 

1 


z — 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


O 


1 

o 


o 


O 


1 

o 


1 

o 


1 

o 


< s: 


CJ 


o 


o 


o 


o 


o 


o 


O 


o 


<-1 


o 


O 


o 


O 


o 


o 


'O 


o 




o 


H- ^ 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


1/1— » 


X 


gj 


>r 


eg 


o 


X 


o 


'J* 


eg 


o 


X 


o 




eg 


o 


X 


sD 


>3- 


eg 


o 


a 


eg 


re 


>r 




o 


gj 




X 


O' 


O 


o 




eg 


re 


sT 


n3“ 


in 


gD 


e- 


OO 




ei 


m 


m 


n 


re 


re 






re 




>r 






>r 


vT 


sT 


>3“ 


nT 


nT 





00 

oO 



u 


U1 


O' 


C-> 


o 


vN 


g> 


X 


O' 


O' 


X 


X 


<T 


n 


n 


O 




X 


o 


re 




1 X 


eg 


o 


O' 




re 


X 


X 


X 


re 


X 


re 


X 


re 




eg 


g5 


o 


n 


u* 


'e 


l/IQ 


O 




«-g 


eg 


re 


re 


■r 


vT 


in 


n 


g) 


X 


e- 




X 


X 


m 


O' 


O' 


o 


z ■« 


in 


in 


in 


in 


m 


in 




m 


in 


in 


in 


in 


in 


in 


n 


n 


n 


n 


in 


>o 


< 








•-g 


«-H 




-g 




—4 


r-H 






, 


—4 






»-H 








CL 


1 


1 


1 


1 


1 


1 


• 1 


1 


f 


1 


1 




1 




1 












h- 


1 

1 


I 

I 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 


1 

1 


1 


1 


1 


1 

1 

1 


1 


1 

1 


1 


1 

1 


\ 

1 


UJ 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

( 


1 

1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


( 


1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

1 


1 

1 


1 

1 


o 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 

1 


1 


1 

1 


1 

I 


1 


1 

1 


z-^ 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


n 


o 


o 


o 


1 

o 


o 


o 


1 

o 


<L ^ 


o 


o 


o 


o 


o 


o 


C_3 


o 


o 


o 


o 


OJ 


o 


o 


o 


o 


CJ 


o 


o 


OJ 




o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


o 


/e 


o 


o 


o 


o 


/e 


o 


o 


1/1 — 


X 


gD 


vt 


eg 


o 


X 


X 


>r 


r\i 




X 


g) 


>r 


r\j 


o 


X 


X 


sT 


eg 


o 




• 


• 




• 




• 


• 




• 


• 




• 




• 














Q 


o 


e- 


X 


O' 


O 


o 


-g 


eg 


re 


<r 




in 


X 


r- 


X 


X 


O' 


o 




eg 


J 










eg 


eg 


eg 


eg 


eg 


eg 


eg 


r\j 


eg 


eg 


eg 


eg 


eg 


re 


re 


*e 



00 

OO 



a 


«-g 




O' 




n 




eg 


in 


re 


O' 




e- 




O' 


re 


•e 


O' 


eg 


eg 


o 


1 ^ 


'4' 

• 


-H 


— < 


eg 


gj 


e- 


n 




g) 


O' 


rg 




>0 


X 


p- 


e- 


gj 


o 


in 


vT 


ncj 


'T 




in 


X 


o 


rsi 


>r 


gj 


r- 


X 


a 




eg 


re 


>r 


in 


g) 


r- 


X 


O' 


Z — 




eg 


C\J 


eg 


-e 


re 


pri 


re 




re 


>3“ 


sT 


nT 


>r 


nT 


vT 


•vT 


sT 




sj- 


< 






•”* 






— 






























a 


1 


1 


1 


1 


1 


t 


1 


1 


1 


1 


f 


1 


1 


1 




1 


1 








H* 


1 

1 


1 

1 


< 

1 


1 

1 


1 

1 


1 

I 


1 

I 




1 

1 


1 


1 

1 


1 

1 


1 


1 

1 


1 

1 


1 

1 


1 


1 


1 


1 

1 


UJ 


1 

1 


1 

1 


\ 

\ 


1 

1 


1 

1 


1 

1 


1 

1 


i 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 


1 

1 

1 


1 


1 

1 


c 


1 


1 


1 


1 


1 


1 


1 


f 


1 


1 


1 


1 


1 


1 


1 

j 


1 


1 

1 


1 

1 


1 


1 


Z-' 


O 


o 


o 


O 


o 


o 


C 


o 


C 


o 


c 


C 


u 


c 


o 


o 


c 


G 


o 


o 


< 2. 


O 


O 


u 


C 


o 


o 


G 


o 


O 


o 


c 


o 


G 


o 


o 


o 


c 


O 


o 


o 


H- 


o 


O 


o 


a 


a 


o 


O 


a 


G 


o 


o 


o 


O 


o 


o 


o 


c 


O 


o 


o 


00 •*“ 


X 


gg 


>3“ 


ev 


o 


X 


>C 


>3“ 


eg 


o 


X 


X 


>3- 


<\ 


n 


X 


sC 


>3“ 


eg 


o 


c 


G 




eg 


re 


'T 


sr 


in 


>C 


e- 


X 


X 


O' 






eg 


eg 


rn 




n 


c 




























— * 

















23 



Figure 4 

Sample TROPOPLOT Output 
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Figure 5 

Sample TROPOPLOT Output 
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transmitter transmission line loss 
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In the process of describing these parameters a brief 

outline of how TROPOPLOT determines path loss is presented. 

A more detailed exposition can be found in Ref. 6. The 

median reference value of attenuation below free space, 

is computed first. The reference value of transmission loss, 

then becomes the sum of the free space attenuation, 

and the reference attenuation, A , 

cr 



L 



cr 



‘-bf 



A 



cr 



dB. 



where the free space loss is defined as. 



= 32.45 + 20 log.|g (f in MHz) + 20 log.|g(d in km). 

The reference attenuation, is determined using one 

of three subroutines depending on the particular mechanism 
of propagation; LOS for line of sight modes, DIFF for 
diffraction, and SCATT for tropo-scatter . Two ray optics is 
used to compute line of sight paths while diffraction paths 
are assumed to be over a double knife edge and a value of 
diffraction attenuation below free space, A^, is computed. 

In the case of scattering the attenuation variable is designa- 
ted as A . The reference attenuation A is determined by 
the smaller value of A^ and A^. In the line of sight case 
the attenuation is calculated for 2 values of distance 
(d^,d.|) for which line of sight propagation is valid to 
produce corresponding values of attenuation A^ and A.j. The 
diffraction attenuation, ‘is computed at the distance 
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and together with and A-j is used to determine the 
slopes, k-| and k.2, of a smooth curve of the reference attenua- 
tion versus distance over the range 1 £ d _< d 

Li b 



Acr " ^0 ^2 ^° 9 ^Q(d/d^) dB. 

The output parameter A^ is defined as 

A = A - k,d - k« log,^d (i.e., d=l ) 
e 0 lo 2 ^lOo' ' 

so for 1 <_d _<d 

Lp 

A = A + k, d + k^ 1 og , rtd 
cr e 1 2 ^10 



When the diffraction attenuation is computed the output 

variables A , and m. come into consideration. The diffraction 
ed d 

attenuation is computed as the weighted average of the 
diffraction attenuation over smooth earth, A^, and the attenu- 
ation over a double knife edge surface, Aj^, where 



A 



j = (l-w)Aj^ + wA^ dB 



A description of the method of weighting these estimators 
and calculating A^ and Aj^ is contained in Annex 3 to Ref. 6. 
The diffraction attenuation, A^, is determined at two 
distances, d^ and d^, in the far diffraction region, and a 
straight line through the points (A^,d^) is 

defined as 



md = (A^ - - *^3) dB/km (slope) 
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and 



^ed " ^fo ^ ^4 ■ '"d‘^4 (intercept) 

where is a "clutter factor" 15 dB). 

If the scatter attenuation is less than the diffraction 

attenuation, as is sometimes the case in trans -hori zon paths 

where the distanced or the angular distance 9 is large, then 

A ^ . When the product of d in kilometers and 9 in 

c r s 

radians is greater than 0.5, A^ is computed at large distances, 
dg and dg, and a straight line through the points (Ag,dg) and 
(Ag,dg) is defined as follows; 

"’s ^ " ^5^ 

and 

^es ^ ^5 ■ "’5‘^5* 

As before the reference attenuation then becomes 



A = A = A^^ 
cr s es 



m d 
s 



d > d 

— X 



The quantity d is defined as the distance where the scatter 
attenuation is equal to the diffraction attenuation. 

It can be seen from the above description and Ref. 6 that 
the actual determination of the output quantities is somewhat 
complex and for most applications the information contained 
in the graphs and tables describing path loss and signal 
strength is sufficient. Consequently the data presented in 
the next sub-section will primarily consist of derivations 
from those portions of the output. If the quantities in 
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Fig. 6 are important to the user it should be noted that for 
certain input combinations the format statements contained 
in the TROPOPLOT program will produce an all asterisk, 
''*******", printout indicating that the format statement 
governing that particular variable requires modification. 

B. ACTUAL PATH LOSS MEASUREMENTS 
1 . Procedure 

In order to obtain the proper input parameters for 
entry into TROPOPLOT certain preliminary measurements and 
calculations were required. The first step was to obtain 
charts of the path area"* and from these determine the terrain 
profile over a great circle path between transmitter and 
receiver. The great circle path was approximated by a rhumb 
line for path lengths less than 70 kilometers and the avail- 
ability of 1 \' topographic charts for the longer paths enabled 
the use of straight lines over the chart area. This profile 
information was plotted on a curved earth's surface using 
the plotting program described in Appendix A. The value of 
the surface refractivity used in determining the effective 
earth's radius was the same as that used in the input to the 



^To obtain a catalogue for U.S. Charts write to 

Denver Distribution Section 
U.S. Geological Survey 
Denver Federal Center, Building 41 
Denver, Colorado 80255 
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program and can readily be obtained from Ref. 9 or from 
actual meteorological data. If the charts of Bean et al . 
(1960) are used this value must be converted from to by 
theequation 



Mg = Nq exp (-0.1057 h^ ) 

The value of h^ used depends on several factors most important 
of which is the mode of propagation. For most of the paths 
considered herein the value of h^ chosen was determined by the 
elevation of the lowest antenna. For trans-horizon paths a 
mean is computed using heights at the obstacle horizons. 

As was previously noted the values of for input can range 
from 250 to 400 N-units, so it is obvious that predictions of 
path loss for anomalous conditions are not readily obtained 
using TROPOPLOT which is designed for long-term and median 
input parameters. The effects of non-standard atmospheric 
conditions along these paths are considered in Section III. 

After plotting the path profile the required input 
parameters were determined using the methods described in 
Appendix B and a number of parameters were varied to produce 
the results described below. 

2 . Results 

The output of the path p rof i 1 e/routi ne for each 
path considered is shown in Figs. 8-19 while the corresponding 
input parameters used are contained in Table I. All the links 
were analyzed for path loss at both VHF and S-band to 
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Figure 8 

for the Zelatched Point to Keyport Path 
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Figure 9 

Terrain Profile for the Makah to Striped Peak Path 
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Figure 10 

Terrain Profile for the Lookout Mt. to Keyport Path 
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Figure 11 

Terrain Profile for the Zelatched Point to 
Bangor Path 
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Figure 12 

Terrain Profile for the Mt. Constitution to 
Gold Mt. Path 
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Terrain Profile for the Striped Peak to 
Mt. Constitution Path 
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Figure 14 

errain Profile for the Bangor to Gold Mt. Path 
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in Profile for the Keyport to Gold Mt. Path 



39 



, 4 



§ 



a 



I 



I 



I 




Vi 






I 

( 









liC* < V 80 AC i£t f€At 




20 40 60 80 ioo 

DISTANCE (km) 



Figure 16 

Terrain Profile for the Lookout Mt. to NOEF, 
Bangor Path 
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Figure 17 

Terrain Profile for the NOEF, Bangor to Gold 

Mt. Path 
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Terrain Profile for the Zelatched Point to JTOEF, 

Bangor Path 
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Figure 19 

Terrain Profile for the Lookout Mt. to Bangor Path 
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illustrate the increase in path loss when frequency is 
increased. The results are shown in Figs. 20-31. The quanti- 
ties (referred to in an above sub-section) which required 
interpretation as to the precise value to be used in detailed 
terrain studies are varied for a single path so that the 
effects of this variation can be shown. A sample path 
(Lookout Mt. to Keyport) was selected for this purpose and 
the effects of several different antenna heights in an 
unmodified version of TROPOPLOT are shown in Fig. 32. These 
same antenna heights are then entered into the modified 
program and the results shown in Fig. 33. To show the effects 
of variations in conductivity and permittivity a and e are 
varied over the same path with the results shown in Fig. 34. 

The effect of variations in the value of Ah is shown in 
Fig. 35. The quantity Ah was computed by taking the mean 
value of Ah over all the paths considered, Ah^ is the tabulated 
estimation of terrain irregularity found in both Ref. 6 and 
Ref. 7, Ah^ is the parameter computed by the methods of 
Appendix B. This comparison is intended to show the degree of 
precision required in choosing a parameter when uncertainty 
over the means of selection is present. 

The importance of detailed path profile information 
is pointed up in the Zelatched Pt. to Keyport path (Fig. 8). 

If TROPOPLOT is used without modification, errors in prediction 

can occur at 2000 MHz the results of which are shown in 

Fig. 20. Note that this output predicts free space attenuation. 
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Figure 20 

Transmission Loss vs. Distance for the Zelatched 
Pt. to Keyport Path 
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Figure 21 

Transmission Loss vs. Distance for the Makah to 
Striped Peak Path 
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Figure 22 

Transmission Loss vs. Distance for the Lookout Mt. 

to Keyport Path 
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Figure 23 

Transmission Loss vs. Distance for the Zelatched 
Point to Bangor Path 
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Figure 24 

Transmission Loss vs. Distance for the Mt. Constitution 

to Gold Mt. Path 
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Figure 25 

Transmission Loss vs* Distance for the Striped Peak to 
Mt. Constitution Path 
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Figure 26 

Transmission Loss vs. Distance for the Bangor to 
Gold Mt. Path 
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Figure 27 

Transmission Loss vs. Distance for the Keyport 
to Gold Mt. Path 
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Transmission Loss vs. Distance for the Lookout Mt. 
to NOEF, Bangor Path 
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Transmission Loss vs. Distance for the NOEF, 

Bangor to Gold Mt. Path 
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Transmission Loss vs. Distance for the Zelatched Point 
to NOEF, Bangor Path 
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Figure 31 

Transmission Loss vs. Distance the Lookout Mt. 
to Bangor Path 
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tenna Height 



Variation for an Unmodified Version of TROPOPLOT 
Transmission Loss vs. Distance 
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Figure 33 

tenna Height Variation for the Modified Version of TROPOPLOT 
Transmission Loss vs. Distance 
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Figure 34 

Variation in z and a for the Lookout to Keyport Path 
Transmission Loss vs. Distance 
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Figure 35 

Variation in Ah for the Lookout Mt. to Keyport Path 
Transmission Loss vs. Distance 
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when there is clearly a large attenuation of the signal by 
path obstacles. The same error can occur when the modified 
program is used; primarily due to the fact that the terrain 
is represented statistically. Large deviations from the 
range of terrain heights occurring in mid-path may not be 
considered as an obstacle by the program when in fact it is 
a major factor in the attenuation process. These smoothing 
effects can be alleviated somewhat by examining the terrain 
profile and performing graphical analysis, such as is 
described in Refs. 1-5, to check the effectiveness of the 
computer routine when its output is questionable. Another 
possible source of error can occur in short paths where only 
a single knife edge obstacle is present. In this case the 
predicted value of path loss is too high since a non-LOS path 
loss is calculated for a double knife edge case. 

3 , Concl us i ons 

It is obvious that the quality of the output of 
this tropospheric propagation prediction program depends on 
which of several possible assumptions were made concerning 
the choice of input parameters. Clearly, if information is 
desired concerning a particular path or path area a path 
profile should be drawn and the ray path checked to insure 
that the values of the corresponding angles of elevation 
(depression) fall within the program constraints. As 
previously stated, if these values exceed the limits then 
TROPOPLOT must be modified and a detailed terrain profile 
obta i ned . 
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The effect of the variations of several of the 
parameters shown above allow some inferences to be drawn 
concerning the care with, which the input variables must be 
chosen. It would seem from the path studied that variations 
of a and e have the least effect on path loss {<1 db), varia- 
tions in Ah have only a small effect, particularly for longer 
distances, and variations in antenna height have the greatest 
effect on path loss. The unmodified version of the program 
appears more susceptible to these changes than the version 
modified to accept detailed profile information. The sensi- 
tivity of the transmission loss to antenna height value 
variation was alluded to in an earlier section and thus the 
suggestion that care be exercised in selecting these values 
is well taken if any degree of precision is desired. The 
detailed computation ofAh, on the other hand, is not necessary. 
By using the tabulated values corresponding to a particular 
type of terrain [Ref. 6] similar results can be obtained. 

The transmission loss predictions seem to be 
accurate, or at least within an expected range of values except 
for short diffraction paths. In order to assess the goodness 
of the prediction, comparison with empirical data for these 
paths should be accomplished in a later study. 
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III. PROPAGATION IN A NON-STANDARD ATMOSPHERE 



The computer analysis conducted in the previous section 
produced a long-term median attenuation value as output and 
consequently required that standard atmospheric conditions 
be assumed. In providing input under this assumption, the 
"standard" atmosphere for the locations under consideration 
was character i zed by the median value of minimum monthly 
surface refractivity at sea level (N^) corrected to emitter 
elevation. This element of the prediction routine could, in 
some cases, lead to an over-optimistic expectation of link 
performance, particularly in a locale where anomalous condi- 
tions frequently occur. One method of avoiding this potential 
source of error is to examine the statistical occurrence of 
these conditions and the severities of the effect of the 
anomalies on the link in question. For the purposes of this 
study the non-standard conditions dealt with involved the 
super-refractive, subrefractive and ducting cases. The 
statistics used in this determination were derived from IREPS 
output and the effects of non-standard refractivity were 
modeled by assuming certain values of effective earth's radius 
for corresponding conditions. Prior to the presentation of 
these results the theory involved in the anomalous propagation 
problem is reviewed. 
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A. THEORY 

The primary medium through which electro-magnetic waves 
with wavelengths less than 1 or 2 meters propagate is the 
troposphere. This region, bounded below by the earth's 
surface and above by the tropopause, is approximately 10 km 
thick and is characterized by a general decrease in tempera- 
ture with height up to a zone of constant temperature called 
the tropopause. The tropopause is not a static boundary, but 
has a height which is variable with both time and latitude. 

The troposphere is usually assumed to be a lossless 
dielectric with y = 1 and a = 0. The index of refraction, n, 
is then 

” *^r atm. 

At the earth's surface n has been found to equal approximately 
1.0003. Since air with a higher water vapor content has a 
larger value of permittivity, that is, 

'^wet air — ^dry air 

and since the lower atmosphere usually has the higher water 
vapor content, the permittivity of the atmosphere exhibits a 
decrease with height to a value of unity. This effect causes 
the refractive index to similarly decrease with height. It 
has been ascertained [Ref. 3] that for temperate climates the 
average variation near the ground is 

4t- = - 0.039x1 0"^ per meter 
dh 
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In order to avoid the use of such small numerical quantities 
a variable, N, designated as the refractivity or co-index of 
refraction, has been defined as 

N = (n-l)xlO^ N-units 

and so within one kilometer of the earth's surface 



_ = _39 N-uni ts/meter 



A standard atmosphere is then defined in Refs. 11 and 12 as 
fol 1 ows : 



1 . £ = 1 ; 
n = 1 



P dT _ 1 °C 

dh 150 m 

77 fi p P = pressure (mb) 

3. N = — (p + 4810 |-) T = absolute temp. (°K) 

e = vapor pressure (mb) 



By assuming a linear refractivity gradient it is possible 
to define an effective earth radius as was done in an earlier 
part of this report. This allows the rays representing the 
radio wave to be drawn as straight lines over a curved earth. 
An alternative repres entati on which enables the effect of 
several different values of effective earth radius to be 
shown on the same profile is to show the rays as curves from 
transmitter to receiver over a plane earth. Details of 
constructing this sort of representation are presented in 
Ref s , 3 and 13. 
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Non-standard propagation can occur when the above condi- 
tions are not met. This is manifested by the manner in which 
the radio wave path is curved. Under normal conditions the 
decrease of the refractive index with height causes a downward 
bending of the wave path. Changes in the refractive index 
super-refraction or in certain cases an upward bending of the 
ray which is termed subrefraction. If the refractive effects 
are severe enough, the formation of ducts occurs in which the 
majority of radio wave energy is trapped within a narrow 
regi on . 

Ducting can occur as a result of a number of different 
atmospheric conditions. The primary requirement in producing 
this phenomena is that the curvature of the radio wave must be 
greater than that of the earth. This is produced by a rapid 
change in N with height which is caused by a sharp decrease in 
moisture with height (abnormal moisture lapse rate) and/or a 
sharp increase in temperature with height (temperature 
inversion). The formation of a ground- or surface-based duct 
can occur when warm dry air flows from land over water. 
Moisture evaporates from the water into the lower layers of 
air cooling the air and producing an increased moisture lapse 
rate and a temperature inversion. The downward motion of warm 
dry air is frequently associated with the clear weather found 
on the eastern side of vast high pressure regions in the lower 
and middle atmosphere. Subsidence such as this can also 
result in elevated ducts which are found to be strongest and 



66 



*» rbf!»' 
•-isqu, 
< Hw V** 

i*2 

'*IOt4(! 
» to 

<r 

■' 3 I 

s-dq *tH| 

9 » •» 1 J t- pH 
QTin?' 

-•JO n»» 

/I 

^ 3 *1 1 



aJII 



lowest in fair weather and highest and weakest near storms. 
Surface ducts are also formed in certain regions by nocturnal 
cooling but unlike the evaporation duct this sort of anomaly 
if found over land, particularly desert regions. 

The opposite effect to super-refractive ducting is the 

I 

upward bending of radio waves caused by an increase in 
refractivity with height. This condition typically occurs 
in moist and cloudy regions. 

i The determining factor in all the above cases is the 

j refractivity gradient. Tab 1 e li be 1 ow , found in Ref. 11, shows 
the range of refractivity gradients over which anomalous 
I condi ti ons prevai 1 . 



Table II. Range of refractive index gradient for differing 
types of propagation conditions. 



p 

(mb ) 


h 

( km) 


Sub- 

Ref racti ve 


Unstratified 


Super- 

Refractive 


1000-850 


0-1 .46 


- d n / d h £ 0 


20£-dn/dh<_60 


100<-dn/dh 


850-700 


1.46-3.01 


-dn/dh<0 


20_<-dn/dh <50 


80^-dn/dh 


700-600 


3.01-4.20 


-dh/dh£0 


20£-dn/dh£40 


70;^-dn /dh 


600-500 


4.20-5. 57 


-dh/dh<0 


20£-dn/dh£30 


50_<-dn/dh 


500-400 


5.57-7.18 


-dn /dh_^0 


20^-dn/dh£25 


40^-dn/dh 



The refractivity gradient can be related to the 
effective earth radius and the resulting ray trace when 
drawn on the terrain profile corresponding to this radius, 
will show the effect of the anomaly on a particular path. 
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The following holds for most conditions [Ref. 3]: 



effective radius = k = 

l+6.4xlO“'^dN/dh 

From this relationship it can be seen that for the standard 
refractivity gradient, which was defined as -39 N-units/km, 
the effective earth radius would be equal to 4/3. The 
presence of anomalous conditions thus cause an effective 
earth radius greater or less than 4/3. The interrelationship 
of the refractivity gradient, the radius of curvature of the 
wavefront, the effective radius, and the assoicated types of 
refraction are shown in Table III [Ref. 3]. 

Note that there is some difference in the exact gradient 

value which defines super-refraction. The super-refractivity 

i 

gradient referred to in Table II includes both the extended 
range and the ducting conditions, while that referred to in 
Table III includes only the ducting case. 

In the actual paths under study several values of 
■ effective earth radius were used to simulate the effect of 
anomalous conditions. Table IV [Ref. 15] provides a guide 

I 

I to the use of k as an estimate of propagation conditions for 
! 99.9%-99.99% path reliability. 
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Table III. Relationship of Refractive Gradient to Various Propagation variaoies. 
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Table IV. K Factor Guide. 









Propagation Conditions 




Perfect 


Ideal 


Average 


Difficult 


Bad 


Weather 


Standard 

atmosphere 


No surface 
layer or 
fog 


Substandard 
light fog 


Surface 
1 ayers , 
ground 
fog 


Fog 

moi s tu re 

over 

water 


Typ i cal 

1 

k 


T emper ate 
zone , no 
fog, no 
ducting, 
good atmos- 
pheric mix 
day and 
night 


Dry, moun- 
tainous , 
no fog 


Flat, 

temperate , 
some fog 


Coastal 


Coas tal , 
water , 
tropical 


F 

i K Factor 


1 .33 


1-1.33 


0.66-1 .0 


0.66-0.5 


0.5-0. 4 



The path profiles from Keyport to Gold Mt., Keyport to 
Lookout Mt., and Bangor to Lookout Mt. are shown for two 
values of K; 0.5 which represents the difficult case as either 
when one antenna is inside a duct and the other outside, or 
the subrefracti ve condition, and -3.57 which represents a 
s ub re f r a c t i ve gradient of -200 N-units/km. The effect of the 
change in earth radius shown in Figs. 36-41 is most noticeable 
in the Keyport to Lookout Mt. path. Using graphical methods 
of analysis found in Refs. 1 and 3 the approximate degradation 
in link performance can be estimated. While little change is 
observed in the Keyport to Gold Mt. link due to the shortness 
of the path length, the Keyport to Lookout Mt. path with K = ^ 
exhibits an additional loss of up to 60 db over that found in 
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Figure 36 

Path Profile for Lookout Mt. to Keyport (K=l/2) 
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Figure 37 

Path Profile for Lookout Mt, to Keyport (K=l) 
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Figure 38 

Path Profile for Keyport to Gold Mt. (K=l/2) 
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Figure 39 

Path Profile for Keyport to Gold Mt. (K=l) 
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Figure 40 

Path Profile for Lookout Mt. to Bangor (K=l/2) 
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Figure 41 

Path Profile Lookout Mt. to Bangor (K=l) 
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standard atmospheric conditions. This same path however 
suffers no apparent degradation for K = 1 as the Fresnel zone 
clearance remains greater than .6. The Lookout Mt. to 
Bangor path degradation was found to be 12-30 db for K = h 
and 6-15 db for K = 1. From these results it can be seen 
that changes in effective earth radius can have a significant 
effect on path loss, particularly as the path length increases 
The profile representation using N data is the most 
commonly used in the literature, however the depiction of 
ducting or trapping types of profiles is often done with 
modified or M-unit profiles where 

/N r a = actual earth radius 

M = N + X 10 

h = height above sea level 

This adjustment to the N-profile is used in the anomalous case 
since straight rays above a curved earth now become curved 
rays above a planar earth and the duct can be more clearly 
shown. The variability of M curves corresponding to a 
particular propagation condition is shown in Table V which 
was derived from Ref. 16. 

The M-curve is thus a transformation in which the 
relative curvature between the normal of the electro-magnetic 
> wavefront and the surface of the earth is unchanged. A duct 
I is formed whenever the M curve has a relative minimum. The 
base of an elevated duct is the height on the M curve below 
the minimum that has the same M value as the relative minimum. 
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When no such height exists, the base of the duct is at the 
ground and a so-called ground based duct exists [Ref. 12]. 
The main disadvantage in the use of M data is that above the 
duct N is grossly o vercorrected . 



Table V. Variability of M curves with type of transmission 
condi ti ons . 



M 




M-Curve 


Type of 
T ransmi ss i on 


Increases 


linearly 


Straight line, with 
positive slope 


Standard 






Slope decreases near 
earth surface with 
rays curving up 


Substandard 






Slope increases near 
surface wi th less 
upward curve 


Superstandard 






Vertical , no curve 


Very great 
coverage 


Decreases 

height 


wi th 


SI ope linear and 
negative rays 
curve down 


Trapping or 
ducting 



^ Having demonstrated the effect of anomalous conditions 
on several propagation paths and considered the theory 
involved the next step is to determine the frequency with 
which these conditions occur. The method chosen to accomplish 
this was to access the IREPS historical file for the geograph- 
ical area in question as described in the following section. 
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B. INTEGRATED REFRACTIVE EFFECTS PREDICTION SYSTEM (IREPS)^ 
IREPS is an experimental system ultimately designed to 
provide the capability for on-board assessment of the effect 
of atmospheric anomalies on sensor performance. The hardware 
portion of the system consists of a mi n i -compu ter with an 
interactive graphics terminal on which the required inputs are 
entered and the results displayed in an easily understood 
graphic format. The information stored in the systems memory 
includes a refractivity library containing long-term meteoro- 
logical statistics on the occurrence of ducting conditions as 
a function of location. This library data is obtained from 
radiosonde stations throughout the world and can be augmented 
by on-scene re f ractome ter or radiosonde data input in order 
to tailor the output to a particular location. In the case 
of the Washington paths, there were only two upper air sounding 
stations in the vicinity of the sites. One is located at 
Tatoosh Island (48°-24'N, 124°-42'W) and the other at 
Quillayute (48°-00'N, 124°-36'W). The historical propagation 
summaries for each location are shown in Figs. 42 and 43. 

Note that while the statistics generated are geared toward 
radar performance, the frequencies to be used in the link are 
included so that some information concerning potential 
anomalies can be derived. The historical meteorological data 
for each station, shown in Figs. 44 and 45, provides informa- 
tion concerning median conditions in the vicinity of the 



^For a detailed description of the system see Ref. 10. 
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sounding station and shows the seasonal and diurnal variation 
in this data. Based on this data, the program then can 
generate a series of profiles. Figures 46 and 47 are M and 
N profiles in the vicinity of the two sounding stations. 

Also available in the program output capabilities is a ray 
trace diagram from a particular emitter. This is useful in 
a qualitative sense in that it shows potential areas of non- 
coverage which may be important in maintaining a ground to air 
communications link. The first set of traces demonstrates the 
effect on the ray trace diagram of selected sites for the 
median ducts found in the Tatoosh Island profile and are 
shown in Figs. 48-52. A second set of traces are shown in 
Figs. 54-57 exhibiting the effect of an elevated duct between 
1100 and 1400 feet. The profile responsible for this condi- 
tion is shown in Fig. 53. The third set of traces show a 
duct between 2300 and 2900 feet in Figs. 59 and 60 and the 
corresponding profile in Fig. 58. 

Besides the above mentioned output format the IREPS 
program can also provide path loss information, with a 
detection threshold for a particular sensor superimposed on 
the distance vs. dB display and also a coverage diagram for 
a particular sensor. These two output diagrams were not 
produced for this study but an example of each is included in 
Fi gs . 61 and 62 [Ref. 10] . 
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Figure 42 

Historical Propagation Conditions for Tatoosh Island, Wash. 
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Figure 43 

Historical Propagation Conditions for QiiilJayute, Wash. 
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Figure 44 

Historical Meteorological Data for Tatoosh Island, Wash. 
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Figure 45 

Historical Meteorological Data for Qulllayute, Wash. 
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LOCATION ‘ 48 24 N 124 36 U 

TINE ALLYEAR DAY 8. NIGHT 

( PLOT , EDIT . LIST , SUNRY , RAYS , LOSS , COOER , END )7 

Figure 48 

Winchelsea Island Ray Trace for Median Conditions 
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LOCATION ^ 48 84 N 184 36 U 

TIME ALL YEAR DAY 8. NIGHT 

( PLOT , EDIT , LIST . SUNRV , RAYS , LOSS , COUER , END )? 

Figure 49 

Lookout Mt. Ray Trace for Median Conditions 
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Figure 50 

Makah Ray Trace for Median Conditions 
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Striped Peak Ray Trace for Median Conditions 
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Figure 52 

Mt. Constitution Ray Trace for Median Conditions 
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( PLOT , EDIT , LIST , SUMRY , RAYS , LOSS , COUER , END )7 



TRf^NSM I TTERy RADAR 
HEIGHT 24SA. 



H 

U 



I 

I 



z 

o 

N 

qe: 

0 

1 

<r 

z 

O 

Z 



> 

(H 

a 

LiJ 

H 

<r 

o 

M 

Q 

Z 




UJ 



o 

M 

f- 

D 

<r 

z 



UJ 

o 

z 

<r 

q: 



XU»-hC!IH t-^z U.UUJH 



Q 

z 

Ui 

01 

Ui 

I> 

o 

o 

(/) 

o 

u 

(0 

> 

a q: 

UI 

M > 

u. q: 

M Q z 
O Ui D 
y M v> 

n y ' 

(/) M H 

o v> 

HU M 

oy J 
z ^ 

• H 

Z O Q 

o z u 



<ru o 

O Z J 

o I-* y 

u H w 



4 J 

a 

3 

Q 

XJ 

(D 

4 J 

cd 

> 

0 ) 

iH 

W 

C 

cd 



o 

M-l 







u 

M-l 

dJ 

a 

cd 

u 

Eh 

cd 



93 




6 . 

r i 

I 







I 




♦ 



STRIPED PK 



l- 

u. 



a 

<r u) 
O U) 
<i OJ 

q: ^ 
\ 

a: • 

UJ • 
f- 

H • 



Z. H 
I 
Z C3 

<r M 

q: Lij 



I 

I 



g I 

N > 

q: 

o Q 
X u 
K 

J <x 

<c o 

Z M 

q: q 

o z 

Z M 




X U t-i O I P 



t-H z U. UJ U H 



Q 

UJ 

M 



u. 

M Q 
O UJ 
UJ M 
Q. U. 
(A t-i 
O 
H UJ 
O 0. 
Z (/> 



- P 

z o 
o z 



p •• 
<x UJ 

o z 

O HI 

J K 



94 



( PLOT . EDIT . LIST . SUNRV , RAYS , LOSS , COUER , END )? 

Figure 55 

Ray Trace from Striped Peak for an Elevated Duct 
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(PLOT /EDIT /LIST/ SUMRV/ RAYS / LOSS / COUER / END )? 

Figure 56 

Ray Trace from Makah for an Elevated Duct 
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( PLOT , EDIT , LIST , SUMRV , RAVS , LOSS , COUER , END )7 

Figure 57 

Ray Trace from Lookout Mt. for an Elevated Duct 
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LOCATION NOT SPECIFIED 
TinE= NOT SPECIFIED 

< PLOT . EDIT , LIST . SUfIRV . RAYS , LOSS . COUER / END )? 

Figure 60 

Ray Trace from Mt. Constitution for an Elevated Duct 
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LOCHTIOM= NOT SPECIFIED 
TinE= NOT SPECIFIED 

( PLOT / ED I T . L I ST , SUNR V . RH VS , LOSS / COUER , END ) ? 

Figure 60 

Ray Trace from Mt. Constitution for an Elevated Duct 
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Figure 61 

Coverage Diagram (IREPS) 
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Figure 62 

Path Loss Display (IREPS) 
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This system is not yet operational and is still under- 
going testing. The value of the output formats in evaluating 
link anomalous propagation potential is clearly obvious. If 
data were available for each of the sites analyzed in the 
study a more exact picture of the statistics of propagation 
would be avail able. 

Based on the available data for the Washington area 
several inferences can be made concerning anomalous propaga- 
tion. For the frequencies in question the probability of 
extended ranges or ducting is between 8 and 10 percent while 
no statistics are presented for be 1 ow normal refraction. For 
the locations from which the soundings are available it can 
be seen that surface ducting occurs most often in the daytime 
during Spring and that the ducts are predominantly less than 
12 meters in thickness. The lowest occurrence of surface 
ducts at Quillayute is in the spring night which also has the 
highest incidence of elevated ducts. At Tatoosh which is a 
more classically marine climate, there are no ducts observed 
to occur during the autumn and winter either at night or 
during the day. The differences in the median refractivity 
at sea level are not significant and both values, 332 for 
Tatoosh and 324 for Quillayute, would produce an effective 
earth radius greater than 4/3, and consequently better than 
average propagation conditions prevail. Until meteorological 
data is available for the actual sites the ray trace diagrams 
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can provide only an indication of the effect of a certain 
profile on the emitter in question rather than actual informa- 
tion concerning real conditions present in the link components. 
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APPENDIX A 



Path Plotting Program 

In order to plot a path profile on a curved earth surface 
either commercially available profile paper must be used or 
the profile can be manually or computer generated. Instruc- 
tions for the construction of path profile paper for various 
values of earth radius can be found in Ref. 13. Computer 
generated profiles can be easily produced if a plotting 
package is available. For the purposes of this study a 
subroutine known as DRAWP was locally available. Other 
alternative plotting methods included the TEKTRONIX 4012 and 
HP 9830 plotting packages, both of which are available at the 
Naval Postgraduate School. 

The program used made use of a transformation which 
allowed the terrain profile to be plotted on a linear graph 
by adjusting the height at a particular distance for the 
effactive curvature of the earth. The profile is drawn by 
the CALCOMP plotter along with the earth surface at sea level. 
The program also fits a line to the terrain data by the method 
of least squares and forms the difference between the value 
of height found on the line and the actual value of the 
terrain height at the same distance. Aside from the actual 
terrain data points the only input parameters necessary to 
obtain this output are either the refractivity gradient or 
the surface refractivity, and the number of points input. 
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The output provided by the program includes a plot of the 
terrain on the effective earth surface, a table of height vs. 
input distance, a list of the terrain heights in kilometers, 
the effective earth radius, and the difference between each 
point on a least square curve and the actual terrain data 
points without the correction for earth curvature applied. 

In order to perform the curve to fit terrain referred to in 
Appendix B, a new set of points must be used taking only 
those the points visible to both transmitter and receiver instead of 
the entire terrain profile. 

The format of the input data deck is presented below: 
First Data Card: 

Column 1: Code "1". surface refractivity will be 

used as i nput 

Code "2". ref racti vi ty gradient will be 
used as input 



Second Data Card: 

Column 1-10: Number of data points 

(right justified integer <_ 200) 

Column 10-20: Refractivity or gradient 
(right justified integer) 



Third Through Last Data Card: 

The distance points in kilometers are listed in order 
using real numbers (i.e., F10.5 format) in columns 
1-10 11-20, .... 61-70. The heights above sea level 

in feet follow directly after the distance listing 
in the same format. A new card is not started unless 
the last distance point was listed in Columns 61-70 

An example of the latter cards is shown below 
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For the table of distance vs. height below 



Distance 

(km) 


Height 

(feet) 


0.0 


420. 


1.0 


500. 


2.0 


10. 


3.0 


125.5 


4.0 


780. 


5.0 


1000. 



The final data cards would be as follows: 



Column: 


1-10 


11-20 


21-30 


31-40 


41-50 


51-60 61-70 


Card 3 


0.0 


1.0 


2.0 


3.0 


4.0 


5.0 420. 


Card 4 


500. 


10. 


125.5 


780. 


1000. 
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APPENDIX B 



INPUT PARAMETER COMPUTATION 

A. DETERMINATION OF DELTA-H 

The program described in Appendix A was used to produce 
from the path geometry the parameters d^-j 2 Ah. 

The computation of Ah is only partially computer-aided by the 
fitting of a straight line to the uncorrected profile heights 
and subtracting the height of the line from the corresponding 
terrain height. This difference is output in the array D(I). 
The interdecile range is then extracted manually from this 
data by finding the value of the lower boundary of the 90th 
percentile and subtracting from it the value of the upper 
boundary of the 10th percentile. An example of the fit of 
this curve to actual terrain data is shown in the main body of 
this study. 

B. DETERMINATION OF HEl ,2 

For line of sight paths the following determination was 
made for effective antenna heights: 

1. A straight line is fitted to the corrected 
terrain profile heights. 

2. The line is defined by 

B(I) = AO + A1 (X(D) 

in the program defined in Appendix A and should 
be modified for earth curvature by 

B'(I) = B(I) - X^/2A 
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3. The points chosen in the actual profile to which 
this line is fitted should be only those points 
which are visible to both transmitter and receiver. 

4. The effective antenna height is then the difference 
between the fitted line and the height of the 
antenna above sea level. 

The procedure described above defines a re fleeting plane 
between the transmitter and receiver. 

If, on inspection of the actual terrain and the least 
square curve, a good fit does not appear to have been obtained 
due to the roughness of the terrain other methods of obtaining 
effective antenna heights are available. If the foreground 
of the antenna represents a good reflecting surface then the 
antenna height over ground can be used, or multiple least 
square curves can be used if the path consists of a number of 
reflecting planes. In the event that terrain characteristics 
defy this sort of analysis, curves and formulae are available 
in Annex III to Ref. 5 which provide alternate methods for 
computing effective height. 

The curve to fit terrain is also used in the determination 
of effective heights for the case of knife edge diffraction 
which is assumed to occur when the path length is less than 
the smooth-earth radio horizon for each antenna. In this case 

two curves are fitted to each of the diffracted paths forming 
a reflecting plane for each antenna. 

For trans-horizon paths the effective heights are obtained 
by considering the actual height of the antenna above sea level 
and subtracting from that height the average height of the 
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terrain between the transmitter and receiver. To illustrate 
this, consider a link with transmitter height above sea level 
h^Q. The heights at N equi -distant points are selected and 
the mean of the central 80% of these values is computed to 
produce h^ where 

.9N 

^t " 78N 2 ^ti 

i = . 1 N 



27 

e.g., N = 31; ^ 23 *^ti ' " O.i.2, ,30 

i = 3 



The effective transmitter height is then 



h 



el 




for h. < h. 

t — to 



If the mean value is greater than the height above sea 
level then the structural height of the antenna above ground 
is used. 

C. DETERMINATION OF DLl ,2 and TEl ,2 

The angular distance 0 is readily obtained from the 

geometry of the profile as shown in Fig. 2. ®02 

are measured and e is calculated in the program as the 

e 

maximum of either e^-j + 0^2 -d^^/a. The distance d^^ is the 

sum of the distances to the obstacle horizons, d^^^ and d |^2 
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D. MODIFICATION TO TROPOPLOT 



With the above means of determining the parameters 

9 1 o» cl,, TROPOPLOT was modified as follows: 
el ,2 ’ LI ,2 ’ 

1. Delete the following statements from the main 
program : 

IF (H16.LE.2.)) GO TO 12 
IF (H1G.GE.2.) .0R.H1G.LE.5.)) GO TO 10 
Z1 = 5.0 
GO TO 11 

10 Z1 = 1 .0+DSIN(3.1415927*HlG/10.0) 

11 HIE = H1G+Z1*DEXP(-2.0*H1G/DH) 

IF (H2G.LE.2.0) GO TO 15 

IF (H2G.GE.2.0.0R.H2G.LE .5.0) GO TO 13 
11 = 5.0 
GO TO 14 

13 12 = 1 .0+DSIN( 3. 141 5927*H2G/10.0) 

14 H2E = H2G+Z2*DEXP(-2.0*H2G/DH) 

15 CONTINUE 

DLl = DLSl *DEXP(-.07*DSQRT(DH/DMAX1 (5.0 ,H1 E) ) ) 

DL2 = DLS2*DEXP(-.07*DSQRT(DH/DMAX1 (5.0,H2E) )) 

TEl = ( .00065/DLSl )*((DLS1/DL1-1 . )*DH-3.077*H1E) 
TE2 = ( .00065/DLS2)*( (DLS2/DL2-1 . )*DH-3.077*H2E) 

2. Insert the following statements 

909 FORMAT (4F10.5) 

READ(5 ,909) DLl , DL2, TEl , TE2 
after the statement 
22 FORMAT (4F10.5). 

Having accomplished these insertions and deletions the input 
card deck as shown in Ref. 7 is changed as shown below for 
user-supplied data only; 
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First Data Card: 



Column 01-10 


"Type of Terrain" - Integer 


Code: 4 


User Terrain Data Used 
No data suppression 


5 


User Terrain Data Used 
No "output parameters" printed 


Column 11-20 


"DB Loss" - Integer 


Code : 0 


No DB data desi red 
(omit 2nd data card) 


1 


DB data desired 


Column 21-30 


"Distance Between Antennas" - Real 
(1-2000) KM 



Second Data Card: 



Column 01-10 


"TX Power Out" (Watts) 
Positive Real Number 


Column 11-20 


"TX Antenna Gain" (db) 
Positive Real Number 


Column 21-30 


"RX Antenna Gain" (db) 
Positive Real Number 


Column 31-40 


"Transmitter Line Loss" ( db ) 
Positive Real Number 


Col umn 41-50 


"Required Receiver Line Loss" (db) 
Positive Real Number 


Column 51-60 


"Receiver Sensitivity" (dbm) 
Positive Real Number 
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Third Data Card: 



Col umn 01-10 


"Surface Re f ra cti vi ty " 
Integer (250-400) N-units 


Column 11-20 


"Surface Conductivity" 
Real (mhos/meter) 


Column 21-30 


"Relative Dielectric Constant" 
Real 


Column 31-40 


"Interdecile Range" 
De 1 ta-H 



Fourth Data Card: 



Col umn 01-10 


"Antenna Polarization" 
Real 



Code: 01.00 Vertical Polarization 
Code: -1.00 Horizontal Polarization 



Column 11-20 


"Frequency" Real 
(20-40,000) 


Column 11-20 


"Calculated Antenna Height of 
Transmitting Antenna" 

Real (0.5-3000.0) Meters 


Column 31-40 


"Calculated Antenna Height 
of Receiving Antenna" 

Real (0.5-3000.0) Meters 



Fi fth Data Card: 



Co 1 umn 01-10 


"Transmitter Obstacle Horizon, DLl " (km) 
Real 


Column 11-20 


"Receiver Obstacle Horizon, DL2" (km) 
Real 


Column 21-30 


"Transmitter Elevation (Depression) 
Angle" (radians) 

Real Positive (Negative) 


Column 31-40 


"Receiver Elevation (Depression) 
Angle TE2" (radians) 

Real Positive (Negative) 
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